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ABSTRACT 

We provide preliminary quantitative evidence that a new solut ion to averaging the observed inho- 
mogeneous structure of matter in the universe (| Wiltshire! l2007al lbT) . may lead to an observationally 
viable cosmology without exotic dark energy. We find parameters which simultaneously satisfy three 
independent tests: the match to the angular scale of the sound horizon detected in the cosmic mi- 
crowave background anisotropy spectrum; the effective comoving baryon acoustic oscillation scale 
detected in galaxy clustering statistics; and type la supernova luminosity distances. Independently 
of the supernova data, concordance is obtaine d for a value of t he Hu bble constant which agrees with 
the measurement of the Hubble Key team of ISandage et al.l ((2006) . Best-fit parameters include a 
global average Hubble constant H Q — 61. 7^ '\ km sec -1 Mpc -1 , a present epoch void volume fraction 
of /vo = 0.76^o9j an< l an a S e 01 the universe of 14. 71^5 billion years as measured by observers in 
galaxies. The mass ratio of non-baryonic dark matter to baryonic matter is S.l^'iL computed with 
a baryon-to-photon ratio that concords with primordial lithium abundances. 

Subject headings: cosmological parameters — cosmology: observations — cosmology: theory — dark 
matter — large-scale structure of universe 



The apparent acceleration in present cosmic expansion 
is usually attributed to a smooth "dark energy", whose 
nature poses a foundational mystery to physics. Our 
standard ACDM cosmology, with a cosmological con- 
stant, A, as dark energy, fits three independent observa- 
tional tests: type la supernovae (Snela) luminosity dis- 
tances; the angular scale of the Doppler peaks in the 
spectrum of cosmic microwave background (CMB) tem- 
perature anisotropics; and the baryon acoustic oscillation 
scale detected in galaxy clustering statistics. In this Let- 
ter we provide preliminary evidence that these same tests 
can all be satisfied in ordinary gener al relativity witho ut 
exotic dark energy, within a model (Wiltshire I2007al lr3) 
which takes a new approach to averaging the observed 
structure of the universe, presently dominated by voids. 

Recently a number of cosmologists have questioned 
whether cosmic acceleration might in fact be an artifact 
of replacing the actual observed structure of the universe 
by a smooth featu reless du s t fluid in Einstein's equations. 
(For a review see iBuchertl ((2007) .) The specifi c solution 
to the averaging problem we investigate here ((Wil tshire 
I2007al) realises cosmic acceleration as an apparent effect 
that arises in the decoupling of bound systems from the 
global expansion of the universe. In particular, gradients 
in the kinetic energy of expansion, and more importantly, 
in the quasilocal energy associated with spatial curvature 
gradients between bound systems and a volume-average 
position in freely expanding space, can manifest them- 
selves in a significant difference in clock rates between the 
two locations. This difference is negligible in the early 
universe when the assumption of homogeneity is valid, 
but becomes important after the transition to void dom- 
inance, making apparent acceleration a phenomenon reg- 
istered by observers in galaxies at relatively late epochs. 
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Galaxies and other objects dense enough to be ob- 
served at cosmological distances are bound systems, 
leading to a selection bias in our sampling of cosmic 
clocks. Since the clock rates within bound s ystems are 
closely tied to a universal finite infinity scale (|EllisHl984l : 
gross variations in cosmic clock rates 
are not directly observable in any observational test yet 
devised. However, relative to observers in bound systems 
an ideal comoving observer within a void would measure 
an older age of the universe, and an isotropic CMB with a 
lower mean temperature and an angular anisotropy scale 
shifted to smaller angles. 

A systematic variation in clock rates between bound 
systems and the volume average, which we will find to 
be 38% at the present epoch, seems implausible given 
our familiarity of large gravitational time dilation effects 
occurring only for extreme density contrasts, such as 
with black holes. However, cosmology presents a circum- 
stance in which conventional intuition based on static 
Newtonian potentials can fail, because spacetime itself 
is dynamical and the definition of gravitational energy 
is extremely subtle. The normalization of clock rates in 
bound systems relative to expanding regions can accu- 
mulate significant differences, given that the entire age 
of the universe has been available for this to occur. 

In this Letter we find best-fit parameters for t he two- 
scale fractal bubble (FB) model ((Wiltshirdl20"07al lb]). The 
two scales represent voids, and the filaments and bub- 
ble walls which surround them, within which clusters of 
galaxies are located. The geometry within finite infin- 
ity regions in the bubble walls is assumed to be spatially 
flat, but the geometry beyond these regions is not spa- 
tially flat. The relationship between the geometry in 
galaxies and the volume-average geometry within our 
present horizon volume is fixed by the assumption that 
the regionally "locally" measured expansion is uniform 
despite variations in spatial curvature and clock rates. 
This provides an implicit resolution of the Sandage-de 
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Vaucouleurs paradox ((Wiltshire! l2007a( ): the "locally" 
measured or "bare" Hubble flow is uniform, but since 
clock rates vary it will appear that voids expand faster 
than walls when referred to any single set of clocks. 

As observers in galaxies, our local average geometry at 
the boundary of a finite infinity region is spatially flat, 
with the metric 



ds% =-dT 2 



\r) [d^ + ^dfi 2 



(1) 



Finite infinity regions are contained within filaments and 
bubble walls. These walls surround voids, where the met- 
ric is not given by |T]) but is negatively curved, with local 
scale factor a v . The average geome try is determin ed by 
a solution of the Buchert equations (Buchcrt 2000) , with 
average scale factor a 3 — / W ia w 3 + / V ia v 3 , where / V i <C 1 
and / w i = 1 — /vi are the respective initial void and wall 
volume fractions at last scattering, when the assumption 
of homogeneity is justified by the evidence of the CMB 
and the Copernican principle. It takes the form 



ds 2 



-dt 2 +d 2 (t) dfj 2 +A(f],t) d(l 2 



(2) 



where th e area function A is defined by a horizon- volume 
average (|Wiltshirell2007af ). The time-parameter t differs 
from the wall-time r of (JTJ) by the mean lapse function 
dt — j(t) dr. The geometry ([2]) does not match the local 
geometry in either the walls or void centres. 

When the geometry dXJ) is related to the average geom- 
etry ([2]) by conformal matching of radial null geodesies 
it may be rewritten 



ds^ = -dr 2 



+ rl(fj,T) dQ 2 



(3) 



where r w = 7 (1 - / v ) 1/3 / W i 1/3 f? w (?7, t). Two sets of 
cosmological parameters are relevant: those relative to 
an ideal observer at the volume-average position in freely 
expanding space using the metric ([2]), and conventional 
dressed parameters using the metric ([3]). The conven- 
tional metric ([3]) arises in our attempt to fit a single 
global metric (JTJ) to the universe with the assumption 
that average spatial curvature and local clock rates every- 
where are identical to our own, which is no longer true. 
One consequence is that the dressed matter density pa- 
rameter, £l M , differs from the bare volume-average den- 
sity parameter, f2 M , according to Q M = 7 3 fi M . 

The conventional dressed Hubble parameter, H, of the 
metric ([3J differs from the bare Hubble parameter, H, of 
((2j) according to 
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Since the bare Hubble parameter characterizes the uni- 
form "locally measured" Hubble flow, its present value 
coincides with the value of the Hubble constant that ob- 
servers in galaxies would obtain for measurements av- 
eraged solely within the plane of an ideal local bubble 
wall, on scales dominated by finite infinity regions. The 
numerical value of H is smaller than the global average, 
H, which includes both voids and bubble walls. Eq. j4|) 
thus also quantifies the apparent variance in the Hub- 
ble flow below the scale of homogeneity. Local measure- 
ments across si ngle voids of the domina nt size, diame- 
ter 30/i _1 Mpc (|Hovle and Vogelevll2004l) . should give a 
Hubble "constant" which exceeds the global average H Q 



by an amount commensurate to H Q — H Q . As voids are 
dominant by volume, an isotropic average will produce 
a Hubble "constant" greater than H Q . This average will 
steadily decrease from its maximum at ~ 30ft. -1 Mpc un- 
til the scale of homogeneit y (~ lOO /i -1 Mpc) is reac hed: a 
"Hubble bubble" feature (|Tomita|[200lHjha et al.ll2007ft . 

We report the results of three indep endent cosmolog - 
ical tests. We use the exact solution ((Wiltshire! l2007bl ) 
to the Buchert equations with boundary conditions at 
the surface of last scattering, z, ~ 1100, consistent 
with observations of the CMB. The luminosity distance, 
d L = 7 _1 a (l + z ) r w, and angular diameter distance, 
d A = d L /(l + z) 2 , are referred to the effective dressed 
geometry ([3j. We take an initial relative velocity dis- 
persion, h r i = 0.99999, between walls and voids, and 
initial void volume fraction, 10~ 5 < / v ; < 10~ 2 , at the 
time of last scattering. The the results are insensitive to 
variations of h r i and / V i for physically reason able priors 
on acc ount of the existence a tracker solution (|Wiltshirel 
l2007bf) to which all solutions tend, to within 1% by red- 
shifts of z ~ 37. The solutions are then effectively speci- 
fied by two independent parameters, which may be taken 
to be the global average Hubble constant, H Q , and the 
present void volume fraction, / v o . 

We have teste d the luminosity distance of the FB 
model against the lRiess et al.l (|2007t ) (Riess07) gold data 
set of Snela and find that for 182 data points and two 
degrees of freedom the best-fit x 2 = 162.7, i.e., a \ 2 of 
approximately 0.9 per degree of freedom, which is a good 
fit. We have performed a Bayesian model comparison of 
the FB model against a flat ACDM model with priors 
55 < H < 75 km sec" 1 Mpc -1 , 0.01 < Q MQ < 0.5. 
This gives a Bayes factor of \nB — 0.27 in favour of 
the FB model, a margin which is "not wo rth more than 
a bar e mention" (Kass and Raftervlll995f ) or "inconclu- 
sive" (|Trottall2007n . Thus the fit of the two models to the 
Riess07 gold data set is statistically indistinguishable. 

The Riess07 gold data set omits data in the "Hubble 
bubble" below redshifts of z < 0.023. In the ACDM 
model, there is no clear theoretical rationale for this; it 
is merely observed empirically that a significant reduc- 
tion in the inferr ed Hubble const ant occurs at the Hub- 
ble bubble scale (|Jha et al.ll2007h . In the FB model the 
Hubble bubble is expected as a feature. 

In Fig. 1 we display the residual difference A/i = 
/i pB — /icmpty, in the standard distance modulus, /i = 
51og 10 (d L ) + 25, of the best-fit FB model from that of 
a coasting Milne universe of the same Hubble constant, 
H Q = 61.7km sec -1 Mpc" 1 , and compare the theoreti- 
cal curve with binned data from the Riess07 gold data 
set. Apparent acceleration occurs for positive residuals 
in the range, z < 0.9. It should be noted that the exact 
range of redshifts corresponding to apparent acceleration 
also depends on the value of the Hubble constant of the 
Milne universe distance modulus used to compute the 
residual. In the FB model the magnitude of the gradi- 
ent of the theoretical residual of Fig. 1 by redshift is less 
than that for comparable ACDM models. This reflects 
the fact that the distance modulus approaches that of a 
Milne universe at late times. 

Statistical confidence limits for the Snela data are dis- 
played as the oval contours in the centre of Fig. 2, in 
the (iJ ,f2 M0 ) parameter space. The dressed density 
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parameter is used here, since it is the one whose nu- 
merical value is likely to be closest to that of a FLRW 
model, and is thus most familiar. Note that ^ M0 — 

±(l-/ v0 )(2 + / v0 ) = i(2 + / v0 ) 3 n M0 (|Wiltshird l2007bl). 

In Fig. 2 we also overplot parameter ranges for which 
two independent cosmological tests have been applied. 
The first test is the effective angular diameter of the 
sound horizon, which very closely correlates with the 
angular scale of the first Doppler peak in the CMB 
anisotropy spectrum. It is often stated that the angu- 
lar position of the first peak is a measure of the spatial 
curvature of the universe. However, this deduction re- 
lies on the assumption that the spatial curvature is the 
same everywhere, appropriate for the FLRW models. In 
the present model there are spatial curvature gradients, 
and we must revisit the calculation from first princi- 
ples. Volume-average negative spatial curvature, which 
accords with tests of ell i pticity in the CMB anisotropies 
(|Gurzadvan et al.ll2005l l2007f ). can nonetheless be con- 
sistent with o ur local observati on of the angular scale of 
the first peak (|Wiltshirdl2007al) . 
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Fig. 1. — The difference in the distance modulus, /i - 
5 log 1 Q(d i ) + 25, with d L in units Mpc, of the FB model with H = 

61.7km sec -1 Mpc -1 , Q M0 = 0.326 from that of an empty coast- 



ing Milne universe, with the same value of H . The IRiess et al.l 
(2007) gold data set of 182 Snela is binned using the criterion 
riiAzi = 5.8, where n is the number of data points, and Azi the 
width of the ith bin. The first bin boundary is set at z = 0.023 as 
"Hubble bubble" points with 2 < 0.023 are exc luded. Our bins ar e 
differ very slightly from those used in Fig. 6 of lRiess et al.l (I2007T ) : 
the single outlier point at z = 1.755, falls in its own bin. This point 
which falls below the theoretical curve is not shown here, but is in- 
cluded in the x 2 anal ysis. We use th e original distance moduli 
reported at http://braeburn.pha.jhu.edu/~ariess/R06/sn_sample 
without the suggested systemati c subtraction ol U.32 mag, as we 
follow the Ccphcid calibration of Sandagc et al. (2006). The boxes 
show the standard statistical errors for the binned data using the 
reported uncertainties, which alrea dy account for luminosity cor- 
rections in the MLCS2k2 reduction l lJha et al.1200711 . The whiskers 
indicate how the residuals move relative to the horizontal axis for 
the 2a limits on H Q with Q M0 = 0.326 fixed: light grey corresponds 
to the 2a upper bound, and dark grey to the 2a lower bound. The 
overlap in these two regions has been coloured black. 



Ideally we should recompute the spectrum of Doppler 
peaks for the FB model. However, this requires consid- 
erable effort, as the standard numerical codes have been 
written solely for FLRW models, and every step has to 
be carefully reconsidered. This task is left for future 
work. The test that we apply here is to ask whether pa- 
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Fig. 2. — la, 2a and 3a confidence limits (oval contours) for 
fits of luminosity dist ances of type la supernovae (Snela) in the 
Riess07 gold dataset (Ricss et al. 2007]) are compared to parame- 
ters within the (Q MQ ,H ) plane which fit the angular scale of the 
sound horizon <5 = 0.01 r ad deduced for WMAP (IBennett et al.l 
120031; ISpereel et ani200l) . to within 2%, 4% and 6% (contours 
running top— left to bottom-right); and to parameters which fit 
the effective comoving baryon acoustic oscillation (BAO) scale of 
104fe~ 1 Mpc observed in ga laxy clustering statistics l lCole et al.l 
[20051: [Elsenstein et al][2005h . to within 2%, 4% and 6% (contours 
running bottom— left to middle-right). 

rameters exist for which the effective angular diameter 
scale of the sound horizon matches the angular scale of 
the sound horizon, 8 = 0.0 1 rad, of the ACDM model, 
as determined by WMAP (jBennett et al.ll2003D . Since 
there is no change to the physics of recombination, but 
just an overall change to the calibration of cosmological 
parameters, this is entirely reasonable. 

In Fig. 2 we plot parameter ranges which match the 
5 = 0.01 rad sound horizon scale to within 2%, 4% and 
6%, using the calculation of the sound horizon given by 
IWiltshirel (|2007al §7.2). The 2% contour would roughly 
correspond to the 2cr limit if the WMAP uncertainties for 
the ACDM model are maintained. As this can only be 
confirmed by detailed computation of the Doppler peaks, 
the additional levels have been chosen cautiously. The 
limits shown have been arrived at assuming a volume- 
average baryon-to-photon ratio in the range i] g ^ = 4.6- 

5.6 x 10~ 10 adopted bv lTvtler et~ai1 (|2000f ) prior to the re- 
lease of WMAP 1. With this range it is possible to achieve 
concordance with lithium abundances, while also better 
fitting helium abundances. This potentia lly resolves an 
anom aly. With the 2003 WMAP1 release (jBennett et all 
l2003h . the baryon-to-photon ratio was increased to the 
very upper range of values that had previously been con- 
sidered, largely due to the consequence for the ratio of 
the heights of the first two Doppler peaks. This ratio 
of peak heights is sensitive to the mass ratio of baryons 
to non-baryonic dark matter - rather than directly to 
the baryon-to-photon ratio - as it depends physically 
on baryon drag in the primordial plasma. The fit to the 
Doppler pea k s requ ired more baryons than the range of 
iTvtler eta l. (2000) admitted, when calibrated with the 
FLRW model. In the FB calibration, on account of the 
difference between the bare and dressed density parame- 
ters, a bare value of Cl BQ ~ 0.03 nonetheless corresponds 
to a conventional dressed value f- B0 ~ 0.08, and an over- 
all mass ratio of baryonic matter to non-baryonic dark 
matter of about 1:3, which is larger than for ACDM. 
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This would certainly indicate sufficient baryon drag to 
accommodate the ratio of the first two peak heights. 

The final set of contours plotted in Fig. 2 relate to the 
independent test of the effective comoving scale of the 
baryon acoustic oscillat ion (BAO), as detected in galaxy 
clust ering statistics (|Cole et al.l 120051 : lEisenstein et al.1 
I2005D . Similarly to the case of the angular scale of the 
sound horizon, given that we do not have the resources to 
analyse the galaxy clustering data directly, we begin here 
with a simple but effective check. In particular, since the 
dressed geometry ([3]) does provide an effective almost- 
FLRW metric adapted to our clocks and rods in spatially 
flat regions, the effective comoving scale in this dressed 
geometry should match the corresponding observed BAO 
scale of 104/i _1 Mpc. We therefore plot parameter values 
which match this scale to within 2%, 4% or 6%. 

The best-fit cosmological parameters, using Snela 
only, are H Q = 61.7+*;^ km sec _1 Mpc _1 and / v o = 

0.76^0 09! with 1(7 uncertainties. The values of the 
mean lapse function, bare density parameter, conven- 
tional dressed density parameter, mass ratio of non- 
baryonic dark matter to baryonic matter, bare Hub- 
ble parameter, effective dressed deceleration parame- 
ter and age of the universe measured in a galaxy are 
respectively: 7 = 1.381±g;g^; ti M0 = 0.125t°;°<>°; 

^MO = 0-33_o.i6) (^MO ~~ ^BoV^BO = ^-2A^ ^0 = 

48.2l^ km sec^Mpc -1 ; q = -0.0428^;^; t q = 

14.7^0 5 Gyr. Statistical uncertainties from the sound 
horizon and BAO tests cannot yet be given, but should 
significantly reduce the bounds on / v q, &> mo etc. 

One striking feature of Fig. 2 is that even if Snela are 
disregarded, the parameters which fit the two indepen- 
dent tests relating to the sound horizon and the BAO 
scale agree with each other, to the accuracy shown, for 
va lues of the Hu b ble co nstant which include the value 
of ISandage et al.l (|2006f h However, they do not agree 
for the values of H Q greater than 70 km sec -1 Mpc -1 
which best-fit the WMAP data ([Bennett et alj 120031 : 
ISpergel et alll2007t ) with the FLRW model. 

The value of the Hubble constant quoted by 

ISandage et al.1 (|2006l ) has been controversial, given the 
14% difference from values which best-fit the WMAP 
data with the A CDM model ([Bennett et al.l 120031 : 
ISpergel et"alll2007| ). However, the WMAP analysis only 
constitutes a direct measurement of CMB temperature 
anisotropies; the determination of cosmological parame- 
ters involves model assumptions. We have removed the 
assumptions of the FLRW model, in an attempt to model 
the universe in terms of the distribution of galaxies that 



we actually observe, with an alternative proposal to av- 
eraging consistent with general relativity. Applied to the 
angular diameter of the sound horizon and the BAO scale 
this leads to different cosmologi cal parameters: ones t hat 
agree with the measurement of ISandage et al.l (|2006l ). 

The combination of best-fit cosmological parameters 
that arises is particularly interesting. The numerical 
value of present void volume fraction, / v o, is identical 
to that of the dark-energy density fractio n, Q An , in the 
ACDM model with WMAP (jSpergel et alJl2007l i. If the 
FB model is closer to the correct description of the ac- 
tual universe, then in trying to fit a FLRW model, we 
appear to be led to parameters in which the cosmologi- 
cal constant is mimicking the effect of voids as far as the 
WMAP normalization to FLRW models is concerned. 
This it does imperfectly, since for a flat ACDM model 
Q M0 = 1 — fi A0 , with the result that the best-fit value of 
fi M0 normalized to the CMB does not match the best-fit 
value of & M0 for Snela with the FLRW model, nor for 
other tests which directly probe Fl M0 - For example, it 
has been recently noted that the values of the normal- 
ization of the primordial spectrum <7 g ~0.76 and mat- 
ter content Q MQ ~ 0.24 implied by WMAP3 are barely 
compatible with the abundances of massive clusters de - 
termined from X-ray measurements (jYepes et alj |2007). 
For the FB model, by contrast, the dressed density pa- 
rameter, fi M0 , includes the range preferred in direct es- 
timations of the conventional matter density parameter. 

The integrated Sachs-Wolfe effect provides a further 
interesting test to be determined. Since the observed 
signal is based on a correlati on to clumped structure 
(|Boughn and Crittenden! l2004h . for large scale averages 
any difference from the ACDM expectation would largely 
depend on the difference in expansion history of the two 
models. However, we might expect foreground voids 
to give anisotropies bel ow the scale of hom ogeneity, for 
which evidence is seen (|Rudnick et al.ll2007| ). 

In this Letter we have offered preliminary quantita- 
tive evidence, via agreement of independent cosmologi- 
cal tests, that the problem of "dark energy" might be 
resolved within general relativity. The differences in 
cosmological parameters inferred in the ACDM and FB 
models - including the average Hubble parameter and 
its variance, the expansion age, dressed matter density, 
baryon-to-photon ratio, baryon-to-dark matter ratio, 
CMB cllipticity - are such that the question as to which 
provides the better concordance model can be answered 
by future observations and new cosmological tests. 
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